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Abstract. We investigate new shapes (in multipole space) of local primordial non- 
Gaussianities in the CMB. Allowing for a primordial isocurvature mode along with the 
main adiabatic one, the angular bispectrum is in general a superposition of six distinct 
shapes: the visual adiabatic term, a purely isocurvature component and four additional 
components that arise from correlations between the adiabatic and isocurvature modes. 
We present a class of early Universe models in which various hierarchies between these 
six components can be obtained, while satisfying the present upper bound on the 
isocurvature fraction in the power spectrum. Remarkably, even with this constraint, 
detectable non-Gaussianity could be produced by isocurvature modes. We finally 
discuss the prospects of detecting these new shapes with the Planck data, including 
polarization. 
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1. Introduction 

The very early Universe being a notoriously difficult territory to explore, it is crucial 
to extract as much information as possible from the measurement of cosmological 
perturbations, via the Cosmic Microwave Background (CMB) and large scale structure 
observations, as a way to test or constrain early universe scenarios. This is why 
primordial non-Gaussianity has been the subject of intense study in the last few years 
(see e.g. [1]). 

So far, WMAP measurements of the CMB anisotropics [2] have set the present limit 
/k>£ai = 32 ± 21 (68% CL) [and -10 < f^l^^ < 74 (95% CL)] on the parameter f^l""^ 
that characterizes the amplitude of the simplest type of non-Gaussianity, namely the 
local shape. Although WMAP data hint at a possible deviation from Gaussianity, one 
must wait for the analysis of the Planck data in order to confirm or infirm this trend. 

A detection of local primordial non-Gaussianity would have a tremendous impact 
on our view of the early Universe. This would indeed rule out all single field infiation 
models, which generate only tiny local non-Gaussianity [3], and give strong support to 
scenarios with additional scalar fields, such as another inflaton, see e.g. [4], a curvaton [5] 
or a modulaton [6], which can easily produce detectable local non-Gaussianity (this type 
of non-Gaussianity is also predicted in other scenarios, for instance the ekpyrotic model, 
see e.g. [7]). 

Interestingly, models with multiple scalar fields open up the possibility to 
generate, in addition to the usual adiabatic fiuctuations, primordial isocurvature 
perturbations, corresponding to fiuctuations in the relative particle number densities 
between cosmological fiuids (whether isocurvature fiuctuations can persist in the post- 
infiation era depends on the details of the thermal history of the Universe). The 
amplitude of such isocurvature modes, which could be correlated with the adiabatic 
one, is now severely constrained by observations of the CMB power spectrum [2]. 

In the present work, we focus on local non-Gaussianity generated by these 
isocurvature fluctuations, previously studied in [8, 9, 10, 11, 12]. Improving on 
these earlier works, we show that the presence of an isocurvature mode, in addition 
to the usual adiabatic one, leads in general to an angular bispectrum that consists 
of the superposition of six elementary components: the well-known purely adiabatic 
bispectrum, a purely isocurvature bispectrum, and four other bispectra that arise from 
the possible correlations between the adiabatic and isocurvature mode. Because these six 
bispectra have different shapes in I -space, their amplitude can in principle be measured 
in the CMB and we have estimated, via a Fisher matrix analysis, what precision on 
these six parameters could be reached with the Planck data, including polarization. 

A natural question is then whether realistic early Universe models could generate 
these new bispectra with detectable amplitude. To answer this question, we consider a 
simple class of models, which generates all six bispectra with various hierarchies between 
their relative amplitudes. Remarkably, some of these models produce detectable non- 
Gaussianity dominated by the isocurvature mode, while satisfying the present upper 
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bound on the isocurvature fraction in the power spectrum. These considerations provide 
a strong motivation to look for these new shapes of non-Gaussianity in the CMB data. 

2. Anguleir bispectra 

Let us thus consider several primordial modes, denoted collectively by X^. We will later 
focus on the case of two primordial modes: the usual adiabatic mode, characterized 
by the total curvature perturbation ^ (on uniform total energy density hypersurfaces) 
and a single CDM (Cold Dark Matter) isocurvature mode S = 3(Cc — Cr), with c and r 
denoting CDM and radiation, respectively. 

The multipole coefficients of the temperature anisotropics (AT/T = 'Yliim^irnYim) 
arc related to the primordial modes via the corresponding transfer functions gi{k), 
so that 

air^^M-iyj 1^ (^Y.^'(^)9i{k?j Y:jk). (1) 
Following [13], we substitute this expression into the angular bispectrum and obtain 

r>mim2m3 = /„ „ „ \ _ r>mim2m3i /qN 
hhh — \"'hmiU'l2m2"'l3m3/ — yiil2l3 '^hhh ^ K'^) 

which is the product of the Gaunt integral 

^r/r"^' = / ^'fi ^1-1 (fi) Yl2m2 (A) ilama (A) (3) 

and of the so-called reduced bispectrum 

Ki2i3 = E (^) V [tlkidki] 9l{h)gi{k,)g^^{k,) 

I,J,K \ / J \i=i ) 

POO 

xB''''{k,,k2,ks) / r^drji,{k,r)ji,{k2r)ji,{ksr) (4) 

'7. 



that depends on the bispectra of the primordial X : 

{X'{k,)X'{k2)X''{ks)) = {27rf5{E,k,)B''''{k^, k^, ks) . (5) 

The reduced bispectrum (4) is here the sum of several contributions, thus generalizing 
the purely adiabatic expression given in [13]. 

In most inflationary models, the "primordial" perturbations X^ (defined during the 
standard radiation era) can be related to the fluctuations of light primordial fields 0", 
generated at Hubble crossing during inflation, so that one can write, up to second order, 

= iV>0» + ^iV„V0»50^ + . . . (6) 

where the d(f)"' can usually be treated as independent quasi-Gaussian fluctuations, i.e. 
{6(j)%k)6(f)\k')) = (27r)3 5"^P5^(A;)5(k + k'), with Ps4,{k) = 27r'^k-%H^/2TTf, a star 
denoting Hubble crossing time. Using Wick's theorem, this implies that the bispectrum 
^iJK -g form 

B'^'^ik,, k2, ks) = X'''''Ps4>{k2)Ps^{k3) + X'''''Ps^{k^)Ps^{ks) 
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+ \''^''P5Ah)Ps^{k2) , (7) 

with X^''^^ = 6°''^6'"^NIj^N^ (the summation over scalar field indices is implicit). Let 
us note that the coefficients A are symmetric under the interchange of the last two 
indices. 

After substitution of (7) into (4), the reduced bispectrum can finally be written as 

I. \^ fI,JKjI,JK 

"hhh - 2^ Jnl "hhh y^) 

I,J,K 

where each contribution is of the form 

POO 

bl;il=3 rHral{r)Pi{r)pf^^{r), (9) 
using (lihh) = [hhh + 5perms]/3!, with 

aiir)^^ J edkji{kr)gl{k) (10) 

/3/(r) = ^1 edkj,{kT)gl{k)P^{k). (11) 

In the /?/, we use the adiabatic power spectrum: = {5''^N^N^)Ps^ = APs^, which 
imphes /i^'l^^ = A^'^^/^2_ ^^^^ 

assume that the coefficients A*"^ are weakly time 
dependent so that the scale dependence of A"^ can be neglected. Note that the elementary 
bispectra (9) depend simply on the power spectrum and the transfer functions whereas 
the main dependence on the early Universe model is embodied by the /nl^- 




Figure 1. Plot of the six elementary bispectra as functions of Z3, for h — I2 = 1300. 



If cosmological perturbations depend on a single scalar field, only the adiabatic 
mode ( exists in (6), and (8) reduces to the familiar adiabatic bispectrum with 
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/^'£^ = (6/5) Jnl = ^lij)/ (N^)"^ . In slow-roll inflation, this coefficient is a combination 
of the usual slow-roll parameters and is thus too small to be detectable. 

However, if perturbations are generated by at least two light scalar fields, an 
isocurvature mode can later coexist with the adiabatic mode. The crucial property, 
then, is that the adiabatic and isocurvature transfer functions, which enter into (10) 
and (11), are very different (as illustrated by the respective Ci = {a^^aim) plotted in 
Fig. 2). Consequently, the angular bispectrum is now the sum of six distinct terms (as 
illustrated in Fig.l), with respective weights /^iS /^f , /^f , /Jf , /Jf , f^f. 

In the particular case where the adiabatic and isocurvature modes depend on two 
disjoint subsets of scalar field fluctuations, the angular bispectrum contains only a purely 
adiabatic contribution and a purely isocurvature one. But, in general, the two modes 
can depend on common scalar field(s), which leads to correlations between the adiabatic 
and isocurvature modes. The four mixed contributions to the angular bispectrum must 
then be taken into account. 



3. Example: a curvaton model 

To illustrate this general situation, we consider a simple class of models based on the 
presence of a spectator light scalar field during inflation, dubbed curvaton [5]. This 
curvaton acquires nearly scale-invariant fluctuations during inflation and, later, behaves 
as a pressureless fluid when it oscillates at the bottom of its potential, before decaying. 

Here, we allow the curvaton a to decay into both radiation and CDM with 
the respective branching ratios 7^ and 7c. Since, in general, CDM can already be 
present before the decay, we define the fraction of CDM created by the decay as 
fc = 7cJlo-/(^c + 7c^<t)) where the Q's represent the relative abundances just before 
the decay. 

As shown in [14], the "primordial" adiabatic and isocurvature perturbations, i.e. 
defined after the curvaton decay, can be written in the form (6), with 

N^^— =— (12) 

n! = hfc - r),N^, = [/,(1 - 2/,) - r] , (13) 
(T* erf 



where r = Sjr^a/ii^: — ^ (1 ~ 1r)^a)] is assumed to be small, since significant 

non-Gaussianities arise only if r <^ 1. 

Let us first discuss linear perturbations. It is useful to introduce the curvaton 
contribution to the total adiabatic power spectrum S = {N^y/[{N^y + {N^y], where 
— H/cf) is associated with the inflaton fluctuation, and = 0. S is directly related 
to the correlation C = P/^^s/ ^/PsPc — sgn(/c — r). The isocurvature-adiabatic ratio, 
given by 
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is constrained by CMB observations [2] to be small (typically a < 0.07) which requires 
either fc — r ot S <^ 1. 

Let us now turn to non-Gaussianities. Using (12), one finds /^l'' = 3S^/(2r). 
This is the dominant contribution in the regime fc — r, the other components being 
suppressed. We thus concentrate on the more interesting case S <C 1 to discuss the size 
of the various components in terms of fc and r, considered as free parameters in our 
phenomenological approach. 

In the regime /c <^ r <^ 1, the purely adiabatic coefficient is the smallest one. The 
other ones are enhanced by powers of (—3) (since N^/N^ — N^^/N^^ — —3): 

/nl = (~3)^/nl^ , /nl'' — > (15) 

where p is the number of "5"' in the triplet {I.J,K}. In particular, the purely 
isocurvature coefficient is enhanced by a factor 27, but with the opposite sign: f^f^ — 
—o?j{2r). All coefficients can be significant if r is sufficiently smaller than cP' . 

In the opposite regime r ^ /c ^ 1, the purely adiabatic coefficient is, once again, 
the smallest one. All the coefficients are now positive and enhanced by factors (3/c/r)P, 
where p is again the number of "5"' indices: 

Note that the enhancement factor is much bigger than in the previous case (15). The 
purely isocurvature coefficient, which dominates, is f-^^ — o? j {2fc) and can be large 
if fc is sufficiently small, while the relative size of the other coefficients depends on the 
ratio r/ fc- 

The above results show that a small isocurvature fraction in the power spectrum 
is compatible with a dominantly isocurvature bispectrum detectable by Planck (e.g. 
a ~ 10^^ and r <^ fc — 10^^ yields f^f'^ ~ 5 x IQ-^). Of course, the relations (15) or 
(16), are specific to the models considered here and would be a priori different in other 
models. It is therefore important to try to measure these six coefficients separately, in 
order to obtain model-independent constraints from observations. 



4. Observational prospects 

To estimate these six parameters, which we now denote /*^*^ , the usual procedure is to 
minimize 

^2 ^ ^^^obs _ ^ ^^obs _ ^ /W^W)), (17) 



mim2m3 
hhh ' 



where B is the angle- averaged bispectrum Bui„i^ = J2m I I 

' \ mi m2 ms I 

the matrix denoting the Wigner-3j symbol. The above scalar product is defined by 
{B,B') = Y.kBi,i2i,B[^i,i,/(Tli,i, with the variance ali^i^ = {Bli^Q - {Bi^^^f ^ 
^khh^h^h^h ^1^12^ = ^ being the number of identical indices among 
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Table 1. Fisher matrix. Only the upper half coefficients are indicated, since the 
matrix is symmetric. 



{^1)^2)^3}- For a real experiment the noise power spectrum is added to the Q in this 
expression. The best estimates for the parameters are thus obtained by solving 

B^^^)p^ = B"''') , (18) 

j 

while the statistical error on the parameters is deduced from the second-order derivatives 
of x^, which define the Fisher matrix, given in our case by Fij = {B^^\ B^^^). 

We have computed this Fisher matrix by extending the numerical code described 
in [16] to include isocurvature modes and E-polarization. This treatment takes into 
account the pure TTT and EEE bispectra, as well as all correlations hke TTE and 
TEE (for more details on the inclusion of polarization see e.g. [17]). We have taken into 
account the noise characteristics of the Planck satellite [15], using only the 100, 143, and 
217 GHz channels, combined in quadrature. Our computation goes up to Zmax = 2500 
and uses the WMAP-only 7- year best-fit cosmological parameters [2]. 

From the Fisher matrix given in Table 1, one finds that the 68 % error on the 
parameters /* is given by 

Af = ^/{F-^)ii = {10, 7, 141, 146, 163, 124} . (19) 

The first two uncertainties are much smaller than the last four. This is due to the severe 
suppression of the isocurvature transfer function at high /, which leads to a saturation of 
the signal to noise ratio for the last four parameters. By contrast, the large / behaviour 
of the first two bispectra is governed by the adiabatic transfer function (less severely 
suppressed), and the error on the first two parameters is thus further reduced at higher 
Z, as shown in Fig. 2. Including polarization has significantly improved (between 2 and 
6 times) the precision on the last four parameters. 

It is also instructive to compare (19) with the naive uncertainties A/* = l/\/F\i — 
{5, 4, 55, 82, 44, 64} obtained by ignoring the correlations, or, equivalently by assuming 
that only one parameter is nonzero. In particular, the contamination of the purely 
adiabatic signal by the other shapes increases the uncertainty, but only by a factor 
2, which is rather moderate (a two-parameter analysis with /*^^^ and f'^^\ assuming 
uncorrelated adiabatic and isocurvature modes, yields errors almost identical to the 
single-parameter ones). Another consequence of the correlations is that an isocurvature 
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Figure 2. The two thick hnes are the adiabatic and isocurvature power spectra, more 
precisely 10^^ x 1(1 + 1)Ci/2tt. The thin hnes represent the uncertainties for the six 
non-Gaussianity coefhcients, as functions of the cut-off Imax- 



non-Gaussianity could be mistaken for a (much smaller) adiabatic one by naively using 
the purely adiabatic estimator: for instance, a purely isocurvature bispectrum would 
give a fake = (Fie/Fn)/^^) ~ IQ-^pl 

5. Conclusion 

In summary, we have introduced a novel analysis of the CMB data, including 
polarization, which relies on the decomposition of a generic local adiabatic-isocurvature 
bispectrum into elementary bispectra. The corresponding amplitudes are independent 
of the details of the early Universe scenario and can be extracted separately from the 
data (the correlations between these parameters and the errors on their measurement 
are contained in the Fisher matrix which we have computed explicitly for the Planck 
experiment). The example of the early Universe scenario discussed here shows that 
a detectable CMB bispectrum dominated by isocurvature modes should be considered 
seriously. Looking for isocurvature non-Gaussianity is complementary to the search for 
an isocurvature component in the power spectrum, which is now standard routine, and 
it would be highly desirable to conduct systematically both types of analysis for the 
future CMB data. 

The present study could be extended in several directions. From the observational 
point of view, it would be useful to investigate the contamination by foreground 
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emissions or secondary effects. It would also be interesting to see how the search for 
isocurvature non-Gaussianities could be improved by combining CMB and large scale 
structure observations, or even taking into account the angular trispectrum [18, 19]. 
From the theoretical point of view, it would be natural to generalize our study to other 
types of isocurvature modes [20] . Moreover, it would be worth undertaking a systematic 
investigation of the possible isocurvature non-Gaussianities predicted by various high 
energy physics scenarios to take advantage of this new observational window. 
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